Carrier Density and Magnetism in Graphene Zigzag Nanoribbons 
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The influence of carrier density on magnetism in a zigzag graphene nanoribbon is studied in a ;r-orbital 
Hubbard-model mean-field approximation. Departures from half-filling alter the magnetism, leading to states 
with charge density variation across the ribbon and parallel spin-alignment on opposite edges. Finite carrier 
densities cause the spin-density near the edges to decrease steadily, leading eventually to the absence of mag- 
netism. At low doping densities the system shows a tendency to multiferroic order in which edge charges and 
spins are simultaneously polarized. 
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Graphene sheets and related carbon based nanomaterials 
have attracted attention recently after seminal experiments 
llll^l revealed novel physics related to their unique electronic 
structurel^l^ngraEhene nanoribbonsfl 11 HJlJi i, E 
[illilllinilHHIBESIilllSlMlH lateral 
confinement leads to size quantization and to one-dimensional 
conduction channels whose properties depend qualitatively on 
edge termination character. Neutral zigzag terminated ribbons 
have attracted particular attention because they have a flat 
band, perflectly flat in simple ;r-band models, pinned to the 
Fermi level. In self-consistent field (SCF) theories, including 
ab initio spin-density-functional theories, the flat band leads 
to robust magnetic order Ferromagnetic alignment of spins at 
the zigzag edges is predicted also in treatments going beyond 
mean field 171 18|]. Although the reliability of SCF theories is 
unceitain and not yet tested experimentally, interest in zigzag 
edge magnetism has remained strong because of potential for 
interesting applications in nano-electronics OlSll . 

Most studies of the electronic structure of zigzag terminated 
graphene ribbons have focused on properties of the neutral 
system or systems with substitutional doping 12511 . We study 
the role of gate voltage induced changes in carrier density, i.e. 
gate doping. A related work in the low caiTier doping regime 
with an additional neutralizing background charge explored 
the possibility of stable non-collinear magnetic states ll26ll . In 
neutral systems, SCF theories predict edge magnetization in 
graphene nanoribbons with opposite spin polarizations on op- 
posite edges ll3.fi [Si- In theoretical studies of locally gated 
zigzag ribbon junctions usually the non-interacting electronic 
structure is assumed |27, 28, 29, 30, 31, 32, 33], neglecting the 
possibility of doping-dependent interaction-driven rearrange- 
ments. In this work we show that gate doping leads to changes 
in charge distribution, spin configuration, and total net spin 
polarization, which are accompanied by important modifica- 
tions in electronic structure. Our study is based on the n- 
orbital Hubbard-model SCF theory for the magnetic proper- 
ties of graphene nanostructures j 1 tI i34i l35i [361 l37ll . in which 
an electron of spin a in site ; experiences a repulsive inter- 
action proportional to the density of opposite-spin electrons 
n,-^. The Hubbard-model SCF theory is broadly consistent 
with DFT calculations when the interaction parameter U is 
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FIG. 1: (Color online) Hubbard model SCF-theory phase diagram 
as a function of doping per length Sn (defined in the text) and rib- 
bon width W for nearest-neighbor hopping yo = 2.6eV, The on-site 
repulsion strength was chosen to have a value U = 2eV which re- 
produce the ribbon band gaps obtained in the LDA-DFT calculations 
in reference |17]. We used 1200 A;-points for Brillouin-zone sam- 
pling. The energetic preference for opposite spins (AF) on opposite 
edges is replaced by a preference for parallel (F) spins at larger dop- 
ing. Above a critical doping Sn ~ 0.7 the SCF calculation does not 
find magnetic states. Solutions at finite doping sometimes (AFb and 
Fb) break the inversion symmetry of the ribbon. When non-collinear 
spin (NC) is allowed canted spin solutions midway between AF and 
F configuration become energetically favored at low doping. 



chosen appropriately. TT-orbital Hartree-Fock theory reduces 
to the Hubbard model when only the on-site Coulomb interac- 
tions are retained. We have chosen to use a Hubbard interac- 
tion parameter U = 2eV which reproduces in the undoped case 
the band-gaps obtained by microscopic density functional the- 
ory in the local density approximation. This value is smaller 
than other estimates ll34ll . but has been adopted with a similar 
motivation in some other recent workll35 . 361. 

The Hubbard model mean-field Hamiltonian for each spin 
a is 
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FIG. 2: (Color online) Total energy differences per edge atom be- 
tween the AF (or AFb) and F, Fb, P states as a function of doping 8n 
for a relatively narrow ribbon with N = S atom pairs per unit cell. For 
low doping AFb type solutions with broken charge symmetry are en- 
ergetically favored over AF solutions although the energy difference 
is very small. The non-collinear (NC) spin solutions are lowest in 
energy in the weakly doped regime. 
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FIG. 3: (Color online) Net spin polarization obtained from the total 
electron spin densities ^ = (n^ — «j) / (n^ for AFb, F and Fb 

solutions as a function of doping. AFb solutions collapse into AF 
solutions with zero net spin polarization for high enough doping. F 
and Fb configurations also progressively lose net spin polarization 
as they approach the non-magnetic P limit. The shaded region rep- 
resents the doping regime where non-collinear solutions are favored 
energetically. 



consist of a nearest neighbor tight-binding term with hop- 
ping 7o = 2.6eV connecting lattice sites / and j, the Hubbard 
term representing electron-electron interactions, and an exter- 
nal potential term accounting for the interaction with the con- 
stant positive background charge proportional to a coefficient 
we choose to be v^xt — —U. Given the uncertainty of predic- 
tions implied by particular versions of SCF-theory, the advan- 
tages of this relatively simple model often outweigh disadvan- 
tages. Because the magnetism in zigzag ribbons is essentially 
one-dimensional, we measure doping 8n in units of the num- 
ber of excess electrons per repeat distance a = 2.46A along the 
edge. The corresponding areal density 5n2D = Sn/W where 
the ribbon width W = \/3Na /2 and is the number of atom 
pairs per ribbon unit cell. 

SCF SOLUTIONS AT FINITE DOPING 

The main players in zigzag edge magnetism are the flat 
band states which occupy one-third of the one-dimensional 
ribbon Brillouin-zone (BZ) and are localized[ll] near the rib- 
bon edges, most strongly so near the BZ boundary \k\ ^ n/a. 
In the undoped SCF ground state, electrons of opposite-spin 
are localized near opposite edges of the ribbon and a gap| 19] 
A oc ly^i separates occupied valence and empty conduction 
band ribbon states. By appealing to particle-hole symme- 
try we can limit our discussion of doping to the «-type case 
in which electrons start filling the conduction band. Dop- 
ing causes charge-density variation across the ribbon and to 
a complicated competition between band and interaction en- 
ergies manifested by the variety of SCF equation solutions 
classified below. We label solutions as AF (opposite) or F 
(parallel) to indicate the relative alignment of spins on oppo- 
site edges. The label NC is used indicate non-collinear spin 
solutions. The label b is applied for solutions which break in- 



version symmetry across the ribbon in a way which will be 
explained in more detail later. Finally we use the letter P to 
designate a paramagnetic state with no local spin-polarization. 
The phase diagram in Fig. ([T]) illustrates the sequence of tran- 
sitions AF NC ^ F ^ Fb ^ P in naiTower ribbons. In 
wider ribbons we find an additional AF state region between 
the F and Fb regimes. 

The total energy per unit cell consist of a sum over all the 
occupied single-particle eigenvalues Ekma labeled with k and 
m the band index divided by A^^: the total number of fe-points 
minus a term to account for the double counting correction in 
the interaction 
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where the occupations n\%"'' are evaluated in the local frame at 
lattice site / where spin is diagonal. Their differences between 
different self -consistent solutions are shown in Fig. (|2]i for a 
particular (A^ — 8) ribbon width when only collinear spin so- 
lutions are considered. In the collinear scheme the energy as- 
sociated with breaking inversion symmetry across the ribbons 
is always small and the main trend is a crossover from antifer- 
romagnetic solutions at small 5n Xo ferromagnetic solutions 
for 5n > 0.04 to non-magnetic solutions for 5n > 0.4. For 
the F type solutions and those with broken charge symmetry 
the system has a nonzero net spin polarization as a function of 
doping density. The doping dependence of spin-polarization 
is illustrated for the same A^ = 8 ribbon width in Fig. (O. 
Each of the solution types identified in Fig. ([T]) is associated 
with particular electronic structure features which are illus- 
trated in Fig. (|4|. For the AF solution, finite doping requires 
that states above the interaction induced gap be occupied. For 
small doping electrons start occupying states near the conduc- 
tion band minima. (See Fig. These additional electrons 



3 













5n = 0.03 




AFb - 







- 






5n = 0.03 

AF - 



7r/2 27r/3 



- 






y 5n = 0.1 

Kf - 



27r/3 



- 








^ (Sn = 0.4 

/ Fb _ 



27r/3 



- 






1 P 



27r/3 



ka 







1 1 1 


1 




1 1 1 


1 




1 1 1 

Til 


1 


0.6 










"i — 






ni 




S 0.5 
















D** — ^ 








Sn = 0.03 






5n = 0.03 






Sn = 0.1 




0.4 




AFb 
1 1 1 






AF 
1 1 1 






F 
1 1 1 





4 8 12 16 

y (A) 



4 8 12 16 

y (A) 



4 8 12 16 
2/ (A) 



1 1 1 1 1 




ni 






Sn = 0.4 






Fb 

1 1 1 






4 8 12 16 

y (A) 



4 8 12 16 

(A) 



FIG. 4: (Color online) Upper row. Band structures corresponding to AFb, AF, F, Fb and P spin collinear solutions of the Hubbard-model SCF 
equations for a zigzag nanoribbon with N = S atom pairs in the unit cell. At finite doping the energy gain due to the gap present in the the AFb 
and AF solutions is reduced, favoring the F solution which does not have a gap. Lower row. Up and down spin electron occupation per lattice 
site in the unit cell across the ribbon. The AFb configuration has broken charge distribution symmetry relative to the ribbon center due to an 
unequal occupation of up and down spin bands. All mean-field bands are invariant under k —k. We show only the portion of the ID BZ 
with states close to the Fermi level. 



suffer a large energy penalty due to the neutral solution band- 
gap and have lower energy when spin-polarized. The resulting 
half metallic solution in the spin-coUinear scheme implies a 
non-zero overall spin polarization in the system and is accom- 
panied by a breaking of charge distribution symmetry around 
the ribbon center This asymmetric charge distribution is a 
combined effect of the net spin polarization and the character 
of the AF solution at the neutrality point, in which electrons 
with opposite spin polarizations are concentrated on opposite 
edges flQ"]. If both up and down spin bands were equally occu- 
pied there would be no charge distribution asymmetry around 
the ribbon center. 

In the low doping regime a non-collinear spin-order that 
continuously bridges the intermediate situation between the 
neutral AF configuration and mostly F configuration at higher 
doping is |26] a possibility. In the version of the Hubbard 
model mean-field theory which allows for non-collinear spin 
denisties we must allow for the possibility that the average 
spin polarization on different lattice sites points in different 
directions |38]. This allows a larger variational space within 
a single Slater determinant approximation and can potentially 
lead to lower energy solutions, but the spin label becomes un- 
defined for each single-particle wave function. We verified 



that non-collinear spin solutions are favored energetically ll26ll 
in the Hubbard model calculations for low doping region and 
that the transition to F configuration happens at doping densi- 
ties typically about 20% higher than when only collinear so- 
lutions are considered. The angle between the spin densities 
on opposite edges and the band structure of the non-collinear 
state are represented in Fig. 5. 

In the intermediate doping regime the total energy is min- 
imized by solutions which are more similar to the F neutral- 
ribbon configuration! 19| which do not have an energy gap, 
and are therefore favored by doping. This transition to F-like 
solutions occurs already at a relatively small value of doping 
5n ~ 0.06. The states that are occupied first at finite dop- 
ing are those near the valley points \k\ — 27z/3a that are[19|] 
spread across the ribbon and therefore control the exchange 
coupling between opposite edges. The W-scaling rules of the 
energy bands near the valley points |19] are consistent with 
the W^^ decay law of the threshold doping at which the tran- 
sition to F-like transition occurs in our numerical phase dia- 
gram. In electronic structures with dominantly F character the 
charge distribution symmetry around the ribbon center is pre- 
served. In this case every occupied states, up or down spin 
and valence or conduction edge band, has a symmetric distri- 
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FIG. 5: (Color online) In the weakly doped region canted spin ori- 
entations develop in order to minimize the total energy when non- 
collinear solutions are allowed. Upper row. Band structure and 
spin resolved electron occupation per lattice for a zigzag ribbon with 
A' = 8 and 8n = 0.05. The occupation and spin polarization at each 
lattice site are represented in a local frame where the spin is diago- 
nal. Lower row. Spin polarization and relative orientation of the spin 
direction between different lattice sites represented with the arrow 
heads. 



bution of electron density around the ribbon center. When the 
doping is sufficiently large, however, we find a broken charge 
symmetry solution that we label as Fb. In this state one of 
the occupied conduction bands has AF (unbalanced across the 
ribbon) rather than F (balanced across the ribbon) character. 
In addition to these solutions, we find that for wide ribbons 
there is an intermediate doping region in which AF solutions 
have lower total energy than the F solutions before the Fb so- 
lution is stabilized. The difference in energy between different 
magnetic solutions is small at intermediate and large doping. 

In the high doping regime the magnetic features of the sys- 
tem progressively disappear as the edge state bands become 
filled. The ribbon is found to turn paramagnetic above a crit- 
ical value that increases with the ribbon width and saturates 
around dric ~ 0.7. Considering that edge localized states in 
the conduction bands with ^-points near 2;r/3a < |^| < n/a 
span approximately 1 /3 of the whole Brillouin zone we find 
that the total amount of doping electrons required to fill com- 
pletely the edge for both up and down spins is 2/3, an amount 
that can be surpassed near the mentioned doping saturation 
limit. 



DISCUSSION 

The AF state of zigzag nanoribbons has the unusual fea- 
ture that inversion symmetry across the ribbon is broken in 
opposite senses in the two spin sub-systems iTt HTqIi . Our cal- 
culation suggest that in low doping regime the system can eas- 
ily develop solutions with a charge density that is distributed 
asymmetrically across the ribbon, creating an interesting and 
unusually strong type of multifeiToic behavior |17, 39] in 
which spin polarization and charge density are coupled. We 
expect that transport properties can correspondingly be ma- 
nipulated in interesting interrelated ways by both external 
magnetic fields and external electric fields directed across the 
ribbon. Edge transport should be strongly suppressed, for ex- 
ample, when a transverse electric field is applied which has 
opposite orientations on opposite ends of a ribbon. 

Above a certain critical doping density, which is inversely 
proportional to the ribbon width we find that the sys- 

tem undergoes a transition to a F configuration in which op- 
posite edges have parallel spin polarizations. When doping is 
increased further the spin-configuration is altered yet again, 
restoring inversion symmetry breaking across the ribbon. In 
this high-doping regime the total magnetic condensation en- 
ergy is small and the energy differences between different 
magnetic configurations is small. Eventually at sufficiently 
high doping the Hubbard-model SCF equations have only 
paramagnetic solutions. 

The Si02 substrates on which exfoliated graphene samples 
are usually prepared have electron density inhomogeneities 
of the order of nfu,c ^ 10 cm ^ that can extend over 
lengths of the order of L ~ Ifim. In the limiting case of rib- 
bons with this same width as the puddle sizes a rough estimate 
of doping per unit lattice constant a within each puddle can be 
evaluated with the product of these two quantities 
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This amount of doping can influence the spin configurations 
in the system and the presence of these random perturbations 
is therefore expected to appreciably weaken the tendency to- 
wards magnetic order, especially for wide ribbons. For this 
reason we should expect a better chances of detecting edge 
magnetism in suspended ribbons which have much weaker 
electron density fluctuations. 

The long-ranged character of the Coulomb interaction, ne- 
glected in the present work, is expected to introduce important 
changes in the details of electronic structure especially in the 
regions in which states with different charge and spin con- 
figurations compete closely. The discrepancies can be more 
acute than in the neutral case because the inadequacy of short 
ranged screening can be more relevant when the occupation 
of each lattice site in the unit cell becomes inhomogeneous as 
we depart from half filling. Nevertheless, it is also likely that 
several qualitative features of the solutions are still correctly 
captured by the Hubbard model and therefore can provide use- 
ful hints on the actual behavior of the magnetic configurations 
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in ribbons as a function of doping. Even though fluctuation 
effects we have neglected may work against the formation 
of long-range order in these ID-magnets, the unusually stiff 
ferromagnetic alignment of the spins predicted by mean field 
theories for zigzag ribbons |16] suggests that magnetic order 
could be possible and should be manifested in some way in 
experiments. 
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